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Abstract 

We calculate the signal-to-noise ratio (SNR) of the stochastic gravitational-wave background in 
an extreme case that its spectrum has a sharp fallofF with its amplitude close to the detection 
threshold. Such a spectral feature is a characteristic imprint of the change in the number of 
relativistic degrees of freedom on the stochastic background generated during inflation in the early 
Universe. We find that, although SNR is maximal with the correct template which is proportional 
to the assumed real spectrum, its sensitivity to the shape of template is fairly weak indicating that 
a simple power-law template is sufficient to detect the signature. 

PACS numbers: 98.80.Cq, 04.30.-w, 04.80.Nn 
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I. INTRODUCTION 



The detection of a stochastic background of primordial gravitational waves (GWs) is an 
exciting challenge. Since the GWs are decoupled from other ingredients of the Universe 
after the Planck time, their detection enables us to probe the early Universe long before the 
recombination, providing much information of the early Universe and high energy physics 
at the epoch when they are generated For example, we can probe the equation of state 
w = p/p in the early Universe using the spectrum of a stochastic GW background 0], and 
lepton asymmetry can be evaluated by investigating effects of neutrino free streaming on 
the power spectrum of primordial GWs [3|. 

One of the most probable sources of such primordial GWs is inflation. Inflation was 
proposed as the most natural solution to the difficulties of the standard big bang cosmology, 
such as the horizon problem and the flatness problem [3]. It also generates primordial 
gravitational waves (tensor modes) [sl as well as the primordial density fluctuations (scalar 
modes) j^. The latter has been observed as the cosmic microwave background anisotropies 
by the cosmic background explorer (COBE) satellite 0], the Wilkinson microwave anisotropy 
probe (WMAP) satellite j^, 0], and so on. On the other hand, the tensor fluctuations have 
not yet been detected and only the upper limit on the ratio of the amplitude of tensor 
fluctuations to scalar fluctuations is obtained as r < 0.55 (95% C.L.) at A; = 0.002Mpc~^ j^. 

The detection of tensor fluctuations (GWs) generated during inflation is very important 
in that their amplitudes can determine the energy scale of inflation directly, while scalar 



fluctuations are sensitive to a combination of the potential energy and its derivative [10 



The energy scale of the inflation str ong ly constrains inflation models proposed so far. Fur- 



thermore, the consistency relation [10| which relates the ratio of the amplitude of tensor 
fluctuations to scalar fluctuations r to the spectral index of tensor fluctuations rij-, if con- 
firmed, would be an extremely important signature of single-field infiation and could help 
to discriminate from other mechanisms for the generation of the spectra such as a curvaton 



mechanism or a cyclic universe [12 



The observational programs for detecting such GWs have been proposed as the next 
generation projects. For example, DECIGO is proposed in Japan [l3| and the big bang 
observer (BBO) at NASA [3]. Since the typical amplitude of a stochastic background of 
primordial GWs is extremely small, the signal to noise ratio (SNR) is not so large even 
if noises are suppressed below the quantum level (ultimate DECIGO). Thus, in order to 
enhance the SNR, search templates as well as cross-correlation analysis are necessary for 



observations of a stochastic GW background 15 1 



Although the spectrum of GWs generated during infiation may well be approximated by 
power-law shape, the present-day spectrum is different from the original one. The amplitude 
of GWs is redshifted by the cosmic expansion and the expansion rate depends on the matter 
content of the Universe. The spectrum is roughly proportional to /~^, /°, / for the modes 
which reenter the horizon in the matter, the radiation, and the kinetic-energy dominated 
phases, respectively, where / is the frequency [l|, 1^. Therefore, templates with (broken) 



power-law shapes are considered and their improvements of SNR have been discussed [17[ . In 
reality, however, the present-day spectrum is never power law nor smooth even if the initial 
spectrum is power law. The effective number of relativistic degrees of freedom changes 



with temperature and these changes leave characteristic features in the spectrum [18|.^ The 



^ The effect of quark gluon plasma phase transition is discussed in [19 1 
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changes of the effective number of degrees of freedom depending on the mass thresholds 
induce relatively rapid changes of the amplitude of the spectrum, rather than the change 
of the power-law index. Therefore, it is important to assess the validity of the use of 
the templates of power-law shape as search templates for a stochastic gravitational-wave 
background. 

In this paper, we consider the frequency ranges of 10"^ — lOHz suitable for proposed future 
detectors of GWs (LISA,DECIGO,BBO). These frequency ranges correspond to the modes 
which reentered into the horizon at the temperature lOOGeV — lO^TeV and the spectrum 
is damped due to the electroweak phase transition. In the next section, we briefly review 
primordial GWs produced during inflation and search templates to detect them. In Sec. Ill, 
we compare SNRs for templates with a simple power-law type to those for templates with 
a rapid change of the amplitude and discuss how effective the former templates are for not 
only detection of such GWs but also probing the change of the amplitude of their spectrum. 
We give discussions and summary in the final section. 



II. DETECTING A STOCHASTIC BACKGROUND OF GRAVITATIONAL 
WAVES 

A. Property of a stochastic background of GWs produced during inflation 

In order to discuss the property of an isotropic and homogeneous stochastic background 
of GWs, one introduces the dimensionless quantity r2g„(/) that is the energy density of 
a GW stored in a logarithmic frequency interval around / divided by the critical energy 
density. In terms of the characteristic amplitude of a GW, he, this quantity is expressed as 

figw(/) = |^/>cr, (1) 

where Hq = 72km/sec/Mpc = 2.3 x 10~^^/sec is the present Hubble parameter. The spec- 
trum of stochastic background of GWs generated quantum mechanically during inflation is 
calculated by the quantum field theory of a massless minimally coupled field in inflationary 
background ^ . However, we obtain some insights into the spectrum shape without detailed 
calculation. Until a mode reenters the Hubble radius, the characteristic amplitude he takes a 
constant value proportional to the Hubble parameter, i/inf, when the mode left the Hubble 
radius during inflation. On the other hand, the amplitude damps as 1/a after the mode 
reenters the Hubble radius. Therefore, the present characteristic amplitude is given by 

Here, to is the present time and tk is the epoch when the mode (with wave number k) 
reentered the Hubble radius, 27r/ = k = a{tk)H{tk). If at tk the Universe is dominated 
by an ingredient with the equation of state w, then (naively) from the Friedmann equation 
we obtain H{tk) oc a(tfc)-^(^+"')/2 and thus / oc a(tfc)"^^+^'"^/^. Therefore the shape of the 
present density parameter VL^^{f) is given by 

f^gw(/)oca(4)i-=^-ocr^(^-3-)/(^+^-). (3) 
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Thus we recover Qg^ oc Z*^, /~^, for w = 1/3, 0, 1, respectively. 

This naive estimate neglects the effect of the change of relativistic degrees of freedom 
during the radiation dominated epoch [l8|. During the radiation dominated epoch, the 
energy density of relativistic particle prad = {''^^ /'iO)g*T^ does not scale as a~^. From the 
entropy conservation 2^ 



—g^sT^o.^ = constant, (4) 

we obtain 

Prad oc g^g's^^a"^. (5) 

Here g^ and (7*5 account for the total number of effective relativistic degrees of freedom. As 
long as the Universe is fully thermalized, these two coincide with each other. Taking account 
of the effective relativistic degrees of freedom, for the mode entering into the horizon during 
the radiation era, we have 

^^gw oc fMtkY oc a{tkYH{tkf oc g,{tk)g.s{tk)-^''' . (6) 
Since g^ and g^s coincide for lOOGeV < T < lO'^TeV, we find Vt^^ oc 



B. Detection method of a stochastic background of GWs 

In general, the GW background signal is expected to be very week and is usually masked 
by the detector noises. To detect such tiny signals, it is practically impossible to detect 
the signal from the single-detector measurement. Thus, we cross correlate the two outputs 
obtained from the different detectors and seek a common signal. We denote the detector 
outputs by Si with 

s,{t) = Kit) + n,{t) , (7) 

where i = 1,2 corresponds to the i-th detector, and hi{t) is the gravitational- wave signal and 
ni{t) is the noise. Then, the cross-correlation signal S is given by multiplying the outputs 
of the two detectors and integrating over the observational time: 



pT/2 .T/2 

S= dt dt' si{t)s2{t')Q{t - 1') , (8) 

J-T/2 J~T/2 



where the filter function Q is introduced to enhance the detectability of the GW signals 
and we take the observation time T to be large enough. We also assume that the statistical 
property of both the GW signal and the noise is stationary, which implies that the argument 
of Q depends only on the time difference t — t' . 

The detectability in the context of stochastic background searches is quantified by the 
signal-to-noise ratio for the cross-correlation signal S 

SNR = ^ (9) 
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Under the assumption that the two different detectors (or output data stream) have no 
correlation of noise in the weak hmit {hi -C Ui), the mean and variance are given by 



o tt2 poo 

(s) = / df\fr'Q^i\f\)imQif), (10) 



207r 

r>00 



OO 



rp poo 

{s')-{src^- dfp,i\f\)P2m\Qif)\\ (11) 

^ J —oo 

where Q{f) is the Fourier transform of the filter function Q{t — t') in the frequency domain. 
Pj(|/|) is the noise power spectrum of i-th detector defined by 

{n,{t)n,{t')) = - / rf/e2-/(*-*')p,(|/|). (12) 

^ J -oo 

7(1/1) is the overlap reduction function which characterizes the reduction in sensitivity to a 
stochastic background arising from the separation time delay and relative orientation of the 
two detectors [2l| . If their orientations are coincident and coaligned without any systematic 
noise correlation between them, the overlap reduction function 7(|/|) becomes constant for 
all frequencies /. When the arms of a detector are separated by 90 degrees, 7(|/|) = 1. 

We would like to determine the functional form of Q{f) which maximizes SNR. In order to 
find such a function, we introduce an inner product (^l-B) for any pair of complex functions 
A{f) and B{f) with weight functions Pi(|/|) ■ Psd/I) m- 

/CO 
dfA*{f)B{f)pm)P2m- (13) 
-oo 

As long as Pi(|/|) is positive for all frequencies, {A\B) satisfies the same properties of 
ordinary inner product in three-dimensional Euclidean space. Using Eqs. fllOl) and (ITT]) and 
the inner product defined in Eq. (fT3|) . the SNR is rewritten as [l5| 



3H^ \ " T 



Vio^v (gig) ' 



7(1/1) ^^gw(|/|) 



\f\'PimP2m 



(14) 



Then, regarding functions Q{f) and 7(|/|) r2gw(|/|)/|/|^Pi(/)P2(|/|) as three-dimensional 
vectors, we find that the filter function maximizing the SNR is given by 

ifi'Piiifmm' ^ ^ 

Here we have neglected an overall normalization constant because the SNR ([9]) is independent 
of it. Note that Q{f) also becomes an even function of a frequency /. Hence, in the following 
we assume that / is positive definite. 

It should be noted that the resultant filter function Q{f) depends on not only known 
functions such as the overlap reduction function and the noise spectrum but also an unknown 
function, that is, the spectrum of GWs flg^{f). Then, it is not until we assume the spectrum 
of GWs that we can determine the maximized SNR and the corresponding filter function, 
which requires us to arrange search templates of the spectrum of GWs. However, since the 
number of templates is limited in practice, such a template does not necessarily coincide 
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with the real spectrum of the GWs background. Thus, it is very important to estimate 
SNRs for such templates when the real spectrum of the GWs background is given. More 
concretely, as explained in the Introduction, the spectrum of GWs generated during inflation 
is in fact neither power law nor smooth. Instead, the amplitude of the spectrum changes 
rapidly due to the changes of the effective number of degrees of freedom depending on the 
mass thresholds. Then, we wonder how much it is justified to use templates with a simple 
power-law type, which have been often considered. If such simple templates have large 
enough SNRs, the number of templates is significantly reduced. On the contrary, in case 
their SNRs are small, we would be able to find a sudden change of the amplitude of the 
spectrum if we use templates which take a sudden change of the amplitude into account, 
albeit the number of the templates becomes much larger. Therefore, we compare the SNRs 
of templates with a simple power-law type to those of templates with a realistic form of the 
spectrum. 



III. EFFECTIVENESS OF SEARCH TEMPLATES 

In this section, we compare SNRs for templates with a simple power-law type to those 
for templates with the change of the amplitude for the LIGO II and the next-generation 
space interferometers such as FP-DECIGO and LISA. Then, we discuss how effective the 
former templates are for not only detection of such GWs but also probing the change of the 
amplitude of their spectrum. 

As the spectrum of GW background, taking into account the change of the effective 
number of relativistic degrees of freedom, we adopt the following simple model in which 
^gw(/) is characterized by the step function:^ 



N 



(16) 



where fc is the critical frequency corresponding to the electroweak phase transition (~ 10 ^ 
Hz) and d is the damping factor that reflects the change of the effective number of relativistic 
degrees of freedom g^,, and is a constant. Since r2gw is proportional to g*^^^ (Eq.(l6])), 
the damping factor due to electroweak phase transition in the standard model of elementary 
particles dgu is dgu — [g*i> TeV)/(7*(lGeV)]~^/^ ^ 0.9. Even in the supersymmetric 
extension of the standard model, only 5'=i,(> TeV) is doubled and the damping factor becomes 
d - 2-^/^dsM ^ 0.7. 

As mentioned in the previous section, the optimal filter is determined by the possible 
spectrum of the stochastic background of GWs. Given the noise power spectrum of the 
detector, we compute the SNRs of the GW spectrum Qg^ using Eq. IHM with 

^^^^ |/|3P,(|/|)P2(|/|)- ^ ^ 

Here, we consider templates with a simple power-law type for the optimal filters 

^^mtcr(/)cxr, (18) 



^ In reality, the change of the spectrum is much milder than the step function 18|. Here we consider the 
extreme case so that the SNR might change most drastically. 
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FIG. 1: The overlap reduction function 7(/) for the Hanford and Livingston, LIGO detector pair. 

where a is a constant. We also consider templates with the change of the amplitude /, that 
is, flmtcrif) = ^gw{f) given in Eq. ( fT6l) for another optimal filter. Then, we compare both 
results and examine whether the template Eq.( |T6l) is effective not only for the detection of 
such GWs but also for probing the change of the amplitude of their spectrum. 

Before going into the detailed computation, we may have a qualitative understanding of 
SNR for a flat Ofiiter- From Eq. (fT7|) for a flat l^fiitor and constant 7(/), a "V"-shaped noise 
density Pi{f) becomes a sharper "A" -shaped optimal filter Q{f) and hence the "bandwidth" 
of Q{f) is very narrow. SNR gets its significant contributions from this bandwidth and does 
not care about the change of the filter beyond that frequency band. So the question is how 
much SNR is improved when fc falls into the bandwidth of Q{f). 



A. Sensitivity of LIGO II to stochastic background of GWs 

First we give the results of SNRs for several templates in the setup of LIGO II. For ground 
based detectors like LIGO II, the sources of noise consist of seismic, thermal, and photon 
shot noises. Then, we use the noise power spectrum of the detecto rs g iving the Figure 1 in 
. The fitting function of this noise power spectrum is given by |23l ] : 



/Max[10-^^(//10Hz)-4 + 10-4^■25(//100Hz)-l•^ 10~^^{f /lO^Ezf] ; 10 < / < 3000Hz, 
00 ; otherwise . 

(19) 

The overlap reduction function ■y{f) is calculated by giving each location for the detector 
pair. Fig. [T] shows the overlap reduction function for the Hanford and Livingston LIGO 
detector pair [l^. Then we compute the SNRs of the GW spectrum fig^ using Eq. ffT^ . We 
take T = 10^ sec and ^In = 10"^''. The results are shown in Figure [2] and Tabled We note 
that SNR is proportional to Qn- 
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FIG. 2: The value (top) and ratio (bottom) of SNRs for rifiitcr = ^gw (solid line in top panel) and 
^filter oc (dotted line in top panel) for LIGO II. We set d = 0.7. 
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TABLE I: SNRs for several templates in the setup of LIGO II. We take T = l(f sec, Qn = W~^, 
and fc = 10"^ Hz. SNR^ax represents SNR for Ofiiter = ^gw with fc = 10~^ Hz in Eq. (fT6|) . SNRflat 
for Ofiitcr oc f^, SNRbpl±i for ^mtcr oc /^^. Note that SNR is proportional to Qn- 

Fig. [2] shows the dependence of the value and ratio of SNR^ax and SNRflat on the critical 
frequency fc with damping factor d = 0.7. Here SNRmax and SNRflat are the value of SNR 
calculated in Eqs. (fT^ and (fT7|) with f^mter = ^gw and flmta oc /°, respectively. We find 
that SNRmax can differ from SNRflat around the range 10 — 50 Hz which is below the typical 
frequency range of the noise spectrum of LIGO II due to factor in Eq. f|T7|) . However, 
the difference is quite small. Comparing with the SNR by a flat spectrum, the improvement 
using the "true" filter Qg^ is at most 2% for d = 0.7. In Tabled we give results for other 
d for /c = 10~^Hz. We find no improvement. The bandwidth of Q{f) is very narrow 
(lOHz < / < lOOHz). SNR gets its significant contributions from this bandwidth and does 
not care about the change of the filter beyond that frequency band. 

We thus conclude that search templates with a power-law form are sufficient in practice 
for the detection of GWs generated during inflation, even though the real spectrum is never 
power law. 
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B. Sensitivity of next-generation space interferometers to stochastic background 
of GWs 



Finally, we also give the results of SNRs for several templates with different forms in the 



setup of next-generation space interferometers such as FP(Fabry-Perot)-DECIGO [2J] and 
LISA 113. 

We consider two detectors forming a starlike constellation. In the pre-conceptual design, 
DECIGO is formed by three drag-free spacecraft, 1000km apart from one another, with 
observation frequency band of around 0.1 — l.OHz. As a result, this starlike configuration is 
identical to the pair of detectors whose arms are separated by 60 degrees with the separation 
2/\/3 X 10^ m in the flat ground. Because the typical wave length of FP-DECIGO is around 
10^ — lO^m, much greater than the separation between detectors, the overlap reduction 
function is almost constant in observation frequency band. Then, we find that 7(/) ~ 0.75 
in the low frequency limit under the starlike configuration in fiat ground. 

The signal processing of FP-DECIGO may adopt the same technique as used in the 
ground detectors. The essential requirement is that the relative displacement between the 
spacecrafts be constant during an observation. Adopting the Fabry-Perot configuration, 
while the arm-length of the detector can be greatly reduced without changing the observed 
frequency range, no flexible combination of time-delayed signal is possible anymore. We 
assume that the output data which is available for data analysis is only one for each set of 
detectors. 

The sources of noise in FP-DECIGO consist of photon shot noise in the photo detector and 
acceleration noise from the drag-free system and radiation pressure noise. Each noise spec- 



trum of FP-DECIGO is given in Photon shot noise is iV^hot = 4.8 x 10-^'^{L/km)-'^{l + 
{f / foY) Hz~^, acceleration noise is A^'acci = 4.0 x 10~'^^(L/km)~^/~'^ Hz~\ radiation pressure 
noise is A'rad = 3.6 x 10~^^/~^[1 -|- (///o)^]~^ Hz~^. Here /o is the characteristic frequency 
given by /o = 1/4JFL with the fineness of JF = 10 and L is the arm-length and we as- 
sume L = 1000km, so that /o = 7.5Hz. The noise spectral density is given by the sum: 
Pi{f) = Ngiiot + A'acci + A^rad- The detailed discussion of noise and instrumental parameters 



of FP-DECIGO is given in (2J]. The noise spectrum of LISA (including white-dwarf binaries 



background [26|) is taken from [27|. The results for FP-DECIGO are given in Fig. [31 Tables 
iniandinil We take T = 10^ sec, Qn = 10"^^ and /o = 7.5 Hz and assume 7(/) = 0.75 as 
mentioned above. The results for LISA are given in Fig. H] and Table IIVI There we take 
T = 10''' sec, ^In = 10~^^, fc = 10~^Hz and assume 7(/) = 0.75 ^ which is maximal for LISA 
since the angle between the two adjacent laser beams in LISA is n/S. 

Fig. [3] (Fig. H]) shows the dependence of the value and ratio of SNRmax and SNRflat on 
the critical frequency fc with damping factor d = 0.7 for FP-DECIGO (LISA). We find that 
both the value and ratio of SNRmax and SNRflat increase around 0.1 — lHz(10~^ — 10~^Hz). 
From Table [ITl we find that the value of SNR in each filter is independent of the damping 
factor d. This is because the contribution of the integration for calculation of SNR Eq. f[T^ 
is the largest around 0.1 — IHz in all frequency ranges. On the others hand. Table UTTl and 
Table HVl show that all SNR is sensitive to the damping factor but is insensitive to the change 
of the spectrum. This is because f2gw is regarded as a fiat spectrum with the damping factor 
d around 0.1 - IHz (lO'^ - lO'^ Hz for LISA) in the case of fc < 0.1 - 1 Hz (lO'^ - 10"^ 



Using more realistic j{f} only reduces the SNR and the bandwidth of Q{f) and does not affect the 
conclusion. 
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FIG. 3: The value (top) and ratio (bottom) of SNRs for f^mtcr = ^ 
^^fiiter oc /° (dotted line in top panel) for FP-DECIGO. We set d 



gw (solid line in top panel) and 
: 0.7. 
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TABLE II: SNRs for several templates in the setup of FP-DECIGO. We take T = 10^ sec, = 
10"^^, and /o = 7.5 Hz. SNR^ax represents SNR for ^mtcv = ^gw with fc = 7.5 Hz (= /o) in 
Eq.m- SNRflat for neater oc f, SNRpl±i for nmcr oc 

Hz for LISA). As a result, although SNRs themselves are enhanced because of much smaller 
noise power spectrum, we find no significant improvement in the use of the filter fig^. 



IV. DISCUSSIONS AND SUMMARY 

Motivated by the fact that the spectrum of gravitational waves has fine structures due 
to the change of the number of relativistic degrees of freedom, we have investigated the 
validity of the use of the templates of power-law shape as search templates for a stochastic 
gravitational-wave background. Comparing the SNR using the template of power-law shape 
and the SNR using the template of step function shape, we find that the resulting SNR is 
insensitive to the change of the amplitude. 

Although we have focused on the change of the spectrum associated with electroweak 
phase transition, our analysis is not limited to it and is easily extended to other frequency 
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damping rate: d 0.9 0.7 0.5 0.1 ^mtcr 
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0.510 ^^mter with a = 

0.417 ^^mter with a = — 1 

0.460 nmtcr with a = +1 



TABLE III: The same as Table [n] except /c = 10^'^Hz in Ea. (fT6]) . 
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FIG. 4: The value (top) and ratio (bottom) of SNRs for rimtcr = ^gw (solid line in top panel) and 
^filter oc (dotted line in top panel) for LISA. We set d = 0.7. 

ranges. 

Our results have both bad news and good news. The bad news is that gravitational 
wave measurements with the amplitude close to their detection threshold are insensitive to 
the detailed structure of the spectrum. Although the spectrum beyond 10~^Hz depends 
on the particle physics beyond 1 TeV, gravitational waves measurements themselves do not 
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TABLE IV: SNRs for several templates for LISA. We take T = 10^ sec, Qn = 10"^^ and fc = W'"^ 
Hz. 



SNRmax 4.59 3.57 2.55 

SNRflat 4.59 3.57 2.55 

SNRpL-i 3.75 2.92 2.08 

SNRpL+i 4.14 3.22 2.30 
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discern the difference: probing SUSY via gravitational- wave observations is not feasible if 
its amplitude is so small that a long time observation is required to achieve sufficient SNR 
as discussed here. 

The spectrum of a stochastic gravitational wave background may not be determined by 
single observations. In order to determine the spectrum (in particular, the spectral index) 
of gravitational waves, multiple observations at different frequencies are required. If such 
observations are realized in the future, the change of the number of relativistic degrees of 
freedom may be detected assuming that the power of the spectrum is known independently. 

The bad news becomes good news at the same time. That is, the templates of simple 
power-law shape are sufficient as search templates for stochastic gravitational waves: no 
detailed transfer function, which is dependent on particle physics models for / > 10~'^Hz, is 
necessary. 
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